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Abstract 

We calculate the conversion rate of high-energy neutrinos propagating in 
constant magnetic field into an electron- W pair {y — > W + e) from the imagi- 
nary part of the neutrino self-energy. Using the exact propagators in constant 
magnetic field, the neutrino self-energy has been calculated to all orders in 
the field within the Weinberg-Salam model. We obtain a compact formula 
in the limit of B <C B cr = m 2 /e. We find that above the process threshold 
fi(th) 2.2-10 16 eV x(B cr /B) this contribution to the absorption of neutrinos 

yields an asymptotic absorption length ~ 1.1 m x(B cr / B) 2 x (10 16 eV/E). 
14.60.Lm, 95.30.Cq 
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I. INTRODUCTION 



The study of creation, propagation, energy loss, and absorption of neutrinos in magnetic 
field is important in several astrophysical contexts and in the early cosmology [1]. Neutrino 
self-energy and dispersion relation is modified in magnetized media [2-4], and processes 
where neutrinos radiate electron-positron pairs [y — > v + e + + e~) [5-7] or gammas [y — > 
v + 7) [8,9,7] have been investigated in the range of energies where it is possible to use an 
effective four-fermion interaction and rates have been obtained in the limits of weak and 
strong magnetic fields. As an example of the importance of macroscopic magnetic field as 
an effective source of energy loss for energetic neutrinos, we recall that the estimate [5] for 
the rate v — > v + e + + e~ in the strong magnetic field near the surface of a neutron star is 
about ten times the rate of pair production in the Coulomb field near the nucleus of metallic 
iron. 

Conversion of neutrinos in WMepton pairs in the presence of magnetic fields [yi — > W + l) 
should be considered when studying the propagation of neutrinos of sufficiently high energies: 
we shall show that in this limit this process gives an important contribution to neutrino 
absorption. A similar process, z/7 — > I W + , has been studied by Seckel [10], who shows 
that at energies above the threshold for W production this process is competitive with v v 
scattering at the same center of mass energies. 

The process we are considering, where an extremely energetic neutrino creates a real W, 
is a second order process in the weak coupling constant g, while the radiation of a e + e~ 
pair through a virtual If or a virtual Z is a fourth order process. Therefore, there is an 
energy above which the conversion into an electron-Pf pair becomes the dominant process. 
In addition, since we use the electro-weak lagrangian and not an effective low-energy theory, 
our result is valid also at very high energies, much higher than the W mass; actually, the 
W-electron decay rate contributes significantly to neutrino absorption only in this limit. 
Notice that eventually the real W decays and that in about 10.5% of cases the final state 
contains a neutrino of the same flavor (e.g., u e — > W + e — > u e e + e"), so that the process can 
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be thought of as the radiation of a lepton pair; in about 21% of cases the final state contains 
a neutrino of different flavor (e.g., v e — > W + e — > v^e), and in the remaining 68.5% of 
cases the final state contains hadrons. 

In this paper we use Schwinger's proper time method [11] to calculate the neutrino self- 
energy in homogeneous magnetic fields and then we obtain the probability of decay into 
^-electron pair by extracting the imaginary part of the self-energy. A similar strategy was 
used by Tsai and Erber [12] to extract the photon pair creation probability from the vacuum 
polarization in intense homogeneous magnetic fields. 

In Section 2 we briefly review the notation and derive the one-loop neutrino self-energy 
in constant magnetic field [2] , in Section 3 we obtain the imaginary part of the self-energy 
and the rate of W^-electron pairs creation in magnetic field. The ensuing discussion and 
conclusions are in Section 4. 

II. NEUTRINO SELF-ENERGY IN A CONSTANT MAGNETIC FIELD 

In this section we review the calculation of the one-loop neutrino self-energy in a homo- 
geneous magnetic field, using the exact fermion and gauge boson propagators in a constant 
magnetic field [2] . We consider a magnetic field with magnitude B pointing along the pos- 
itive z-direction. The only non-vanishing components of the electromagnetic field strength 
tensor F^ v are F 12 = —F 21 = B. The exact expression for the electron S(x',x") [11,13] 
and W boson G fJ,u (x',x") [2] propagators in a constant magnetic field are obtained using 
Schwinger's proper time method: 




(1) 





where the translationally invariant parts of the propagators are 
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In the rest of the paper we shall drop the infinitesimal imaginary contribution to the masses 
— ie, which determines the correct boundary conditions; if necessary it can be easily reintro- 
duced: m 2 —> m 2 — ie and M 2 — > M 2 — ie. In our notation, — e and m are the charge and 
mass of the electron, M the VT-mass, x the gauge parameter, a 3 = a 12 = |[7 1 ,7 2 ], and the 
metric is = diag(— 1, +1, +1, +1). We choose the electromagnetic vector potential to be 
= —7;F fMU x l/ and, therefore, the phase factor in Eqs. (1) and (2), which is independent of 
the integration path, is [13] 
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For any 4-vector a M we use the notation = (a , 0,0, a 3 ) and a,j_ = (0,a 1 ,a 2 ,0); it is 
easy to show that 
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Note that the W and Goldstone scalar propagator were obtained in Ref. [2] by introducing 
a new gauge fixing term (EGF) which is manifestly invariant under electromagnetic gauge 
transformations. The advantage of the EGF gauge is that the electromagnetic potential has 
no cross couplings with the W and Goldstone fields and, therefore, these two fields do not 
mix in the presence of a magnetic field. 

In the remainder of this paper we focus our attention on electron-type neutrinos, the 
generalization to // and r-neutrino is straightforward. For the purpose of this work, it 
would seem convenient to work in the unitary gauge (x — > oo), where the unphysical scalars 



disappear. However, the PF-propagator is quite cumbersome in this gauge. We prefer to 
work in the Feynman gauge (x — 1), where the expression of the propagator is much simpler. 
In principle the choice of the Feynman gauge carries the price of calculating an additional 
bubble diagram: the one with the Goldstone scalar. But this scalar bubble diagram does not 
contribute to leading order, since it is suppressed by a factor of m 2 /M 2 xs 4.04 • 10~ n , and 
can be neglected. Therefore, we only need to calculate the bubble diagram with a W-boson: 

Z w (x', x") = l -^- lR ^S(x', x") lvlL G^{x\ x") (7) 

where 

1 + 75 1-75 , fl x 

This expression is translationally invariant, since the W and the electron carry opposite 
charge and therefore the phase factor in the W-propagator cancels the phase factor <fi* in 
the electron propagator. We can write Y^wip) in momentum space using 

Z w (x',x") = J^e^'-^ip) (9) 



as 
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q = p — k , Z\ = eBs\ , z 2 = eBs 2 ■ (11) 

We now do the straightforward 7- algebra, change variables from to z iy translate the k 
variables of integration as follows 

(*« • *J - (h + ■ k ± + tall ^ D ^ ) (12) 

and, finally, perform the four gaussian integrals over the shifted variables k. The result is: 
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and the appropriate counter-terms (c.t.) are defined such that 
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The 5-independent counter-terms are unimportant for the purpose of this work, since they 
do not contribute to the imaginary part of the self-energy 

In order to more easily isolate the range of variables that give most of the contribution 
to the integrals and the terms in the integrand that are negligible, we also find it convenient 
to change integration variables from (-21,-22) to (z,u) 

M 2 z 2 1 -2 2 
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where we have introduced the two parameters (3 = eB/m 2 and rj = m 2 /M 2 = 4.0376 • 10 u ; 
the resulting expression is: 
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III. RATE OF IT-ELECTRON PAIR CREATION 



Because of the oscillating phase exp (—iz), the main contribution to the integral over z 
comes from the region where -2^1. If we are only interested in neutrinos that travel through 
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a "moderate" magnetic field eB <C m 2 = eB cr , this means (5 <C 1, we can expand the terms 
in the integrand in power series of f3z <C 1. Furthermore, since T, w is quite small, of order 
g 2 eB/M 2 , as it can be inferred from the expression (17) or from the explicit calculation in 
Ref. [2], we can set $ = in T, w . 
After this is all done, we obtain 

^2 /„E?\ 2 ™ rM 2 Im? r2 \ 



g feB\ Z" 00 rM 2 /m 2 
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Notice that the term in the exponential proportional to z 3 cannot be dropped in general, 
in spite of the fact that (3z <C 1: the coefficient of (f3z) 2 , i.e. (mp±/M 2 ) 2 , could be very large 
when the neutrino is very energetic. In fact, there are two physically interesting regimes 
that are discriminated by the dimensionless field dynamical parameter 

• < 2 °> 

which can be read from the ratio of the z and z 3 term in the exponential of Eq. (19) or 
inferred from kinematical considerations and the momentum change of a virtual e — W pair 
in a magnetic field. 

At low energies, ^ C 1, the z 3 term in the exponential can be dropped, the self-energy 
is real and we obtain the result of Ref. [2]; at high energies, £ >> 1, the self-energy acquires 
an imaginary part that we shall calculate in the following. In this last case we write Tiw as 

^ - Iwy i§) ^/>/ >--<* + » + ^ v ' (21) 

where we have dropped all nonleading terms in rju and extended the integration in du to 
oo, due to the fact that the main contribution to the integral in du comes from the region 
u ;$ 1 because of the oscillating phase exp (—izu) and that rj = m 2 /M 2 is extremely small. 

The integral in dz of the imaginary part of the self-energy can be performed with the 
substitution z = yy/1 + uj (£u) in terms of the modified Bessel function 
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obtaining 
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The final integration in du yields 

and, therefore, the absorption coefficient, a = — 2p±QT l w(p), is 

a ~ UKhchcKM) Kb J [ +v m P± b) 6 • [Zb) 

From the exponential in Eq. (25) we can read the threshold for the process 1 : if we define 
= y/s then 

CJ9 f) = v ^Mc 2 ^^^2.2-10 16 eV^) . (26) 
For energies well above the threshold, the absorption coefficient has the asymptotic behavior 
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where we have used the numerical value G F /(hcf = g 2 /[V2(2Mc 2 ) 2 hc] = 1.16639 • 
l(T 5 GeV- 2 . 

In case we consider (z/ T ) instead of v e the corresponding threshold energies are higher 
by a factor m^/m e ps 206.8 (m T /m e ps 3478), while the asymptotic behavior does not change, 
since it depends only on M and not on m 2 . 



1 The kinematical threshold is obviously given by the sum of the masses M + m, below which the 
rate is rigorously zero: this threshold is lost in our expansion, but it is unimportant as long as the 
effective threshold, below which the process is exponentially suppressed, is very much larger than 
M + m. 

2 In principle the full result would have additional dependences on the lepton mass coming from 
terms that contain r\ = (m/M) 2 and that we have disregarded being extremely small for electrons; 
these corrections are larger for muons and especially taus, but still small. 
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IV. DISCUSSION AND CONCLUSIONS 



In Figure 1 we show the absorption coefficient for v e — > W + e as function of the neutrino 
transverse energy E = p_i_c for several values of the magnetic field B: B = 10~ 1 B cr , B = 
10- 2 B cr , B = 10- 3 B cr , B = 10- A B cr , and B = 10- 5 B cr with B cr = m 2 /e = 4.4 • 10 10 Gauss. 
It is evident that the process has an energy threshold that grows for smaller values of B, 
as quantitatively described by Eq. (26). If we consider the angular coefficient, which is one, 
and the spacing between the curves above the threshold, we find that a grows linearly with 
E and quadratically with B, in agreement with Eq. (27). The same Figure 1 is valid for 
(v T ), if one multiplies both the horizonal and the vertical scale times (m M /m e ) = 206.768266 
((m T /m e ) = 3477.6). 

In Figure 2 we compare the process u e — > W + e with the process v — > v + e + + e~ that 
has a much lower threshold; this threshold can be estimated from the dimensionless field 
dynamical parameter characteristic to this process k = eBp±/m 3 , which is analogous to £ 
of Eq. (20) with the substitution M — > m: the threshold is smaller by a factor of about 
m 2 /M 2 ~ 4 • 10 -11 . The curves plotted as function of k show clearly that the thresholds 
of both processes are function only of the product Bp±, see Eq. (26). For the process 
v — > v + e + + e~ we use the result shown in Eq. (8) of Ref. [6]: 

+ x G 2 f(9v + 9a) B\ 2 (, . . ln(3) 29\ . . 

a(u -> ue + e~) = ^—^ m 4 E — ln(«) - -ii - 7e - - , (28) 



(3tt) 3 V5 cr 7 V 2 24 / 

which the authors claim valid for E < M 3 /eB (k < (M/m) 3 w 4 • 10 15 ): we plot Eq. (8) 
of Ref. [6], see Eq. (28), up to k — 10 15 . We show results only for B = 0.1B cr (top curves) 
and B = 0.001B cr (bottom curves). We see that the process v — > v + e + + e~ (dashed 
curves) dominates below the threshold of z/ — > W + e (solid curves), but above this threshold 
z/ — > + e is almost two orders of magnitude larger (about a factor 50). Therefore, above 
the threshold E pz (mM 2 )/eB (k (M/m) 2 m 2 • 10 10 ), the effective lagrangian (four 
fermion interaction) cannot be used. The rate of v — > z/ + e + + e~ can instead be estimated 
using the rate of v — > + e, which gives the total adsorption, times the branching ratio of 
W -> i/ e + e, which is (10.66 ± 0.20)%. 
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In conclusion, we have calculated the absorption rate of very high energy neutrinos in 
magnetic field. Our main result is given by the compact formula in Eq. (25) valid for electron 
neutrinos. The result for muon or tau neutrinos can be obtained by substituting the electron 
mass m with the muon or tau mass remembering that the electron mass m appears also in 
B cr = m 2 /e. 

This process is exponentially suppressed, and, therefore, it can be disregarded, for en- 
ergies below a threshold energy inversely proportional to the magnetic field; for a field one 
tenth of the critical field this energy is of the order of 10 17 eV, see Eq. (26) and Fig. I. 

Above this threshold the absorption coefficient grows linearly with energy and quadrat- 
ically with the field as shown in Eq. (27): for a field one tenth of the critical field and an 
energy of 10 18 eV the absorption coefficient is about one inverse meter (see Fig. 1). 

Above the threshold this process substitutes the radiation of e + e~ pairs as dominant 
mechanism for v absorption in magnetic field (see Fig. 2). 
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E [eV] 

FIG. 1. Neutrino absorption coefficient for the process v e — ► W + e in inverse meters as func- 
tion of the neutrino transverse energy (E = p±c) in eV for five values of the magnetic field: 
B = lCT 1 ^, B = 10~ 2 B cr , B = 10~ 3 .B cr , B = W~ 4 B cr , B = lQ~ b B cr going from top to bottom, 
with B cr = 4.4 • 10 10 Gauss. The same curves apply for iy T ), if both the horizonal and the 
vertical scale are multiplied times (m M /m e ) = 206.768266 ((m T /m e ) = 3477.6). 



12 




FIG. 2. Neutrino absorption coefficients for the processes v e — > W + e, (solid curves) and 
z/ — > f + e + + e~ (dashed curves) in inverse meters as function of the dimensionless characteristic 
parameter k = eBp±/m 3 for two values of the magnetic field: B = W~ 1 B cr and B = 10~ 3 B cr 
going from top to bottom. 
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